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Abstract: We present a detailed analysis of the symmetry properties of a four-quark wave 
function and its solution by means of a variational approach for simple Hamiltonians. We 
discuss several examples in the light and heavy-light meson sector. 
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X 

e3 ' 1. Introduction 

_ _ i 

The potentiality of the quark model for hadron physics in the low-energy regime became first manifest 
when it was used to classify the known hadron states. Describing hadrons as qq or qqq configurations, 
their quantum numbers were correctly explained. This assignment was based on the comment by 
Gell-Mann [1] introducing the notion of quark: "It is assuming that the lowest baryon configuration 
(qqq) gives just the representations 1, 8 and 10, that have been observed, while the lowest meson 
configuration (qq) similarly gives just 1 and 8". Since then, it has been assumed that these are the 
only two configurations involved in the description of physical hadrons. However, color confinement is 
also compatible with other multiquark structures like the tetraquark qqqq first introduced by Jaffe [2] . 
During the last two decades there appeared a number of experimental data that are hardly accommodated 
in the traditional scheme defined by Gell-Mann. 
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One of the first scenarios where the existence of bound multiquarks was proposed was a system 
composed of two light quarks and two heavy antiquarks (nnQQ). These objects are called heavy-light 
tetraquarks due to the similarity of their structure with the heavy-light mesons (nQ). Although they may 
be experimentally difficult to produce and also to detect [3] it has been argued that for sufficiently large 
heavy quark mass the tetraquark should be bound [4, 5]. The stability of a heavy-light tetraquark relies 
on the heavy quark mass. The heavier the quark the more effective the short-range Coulomb attraction 
to generate binding, in such a way that it could play a decisive role to bind the system. Moreover the QQ 
pair brings a small kinetic energy into the system contributing to stabilize it. 

Another interesting scenario where tetraquarks may be present corresponds to the scalar mesons, 
J PC = ++ . To obtain a positive parity state from a qq pair one needs at least one unit of orbital angular 
momentum. Apparently this costs an energy around 0.5 GeV 1 , making the lightest theoretical scalar 
states to be around 1.3 GeV, far from their experimental error bars. However, a qqqq state can couple to 
J PC = ++ without orbital excitation and, as a consequence, they could coexist and mix with qq states 
in this energy region. Furthermore, the color and spin dependent interaction arising from the one-gluon 
exchange, favors states where quarks and antiquarks are separately antisymmetric in flavor. Thus, the 
energetically favored flavor configuration for qqqq is [(qq)3(qq)3\, a flavor nonet, having the lightest 
multiplet spin 0. The most striking feature of a scalar qqqq nonet in comparison with a qq nonet is a 
reversed mass spectrum (see Figure 1). One can see a degenerate isosinglet and isotriplet at the top of 
the multiplet, an isosinglet at the bottom, and a strange isodoublet in between. The resemblance to the 
experimental structure of the light scalar mesons is striking. 



Figure 1. Quark content of a qq nonet (left) and a qqqq nonet (right). 
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Four-quark states could also play an important role in the charm sector. Since 2003 there have 
been discovered several open-charm mesons: the _D*j(2317), the D s j (2460), and the -Dq(2308). In 
the subsequent years several new states joined this exclusive group either in the open-charm sector: the 

1 This effect can be estimated from the experimental M(L = 1)—M(L = 0) energy differences: ai(1260)— p(776) = 484 
MeV, /i(1282) - w(782) = 500 MeV, hi (1170) - ?7(548) = 622 MeV, h c (3526) - r/ c (2980) = 546 MeV, Xd(3511) - 
J/*(3097) = 414 MeV, Xbi(9893) - T(9460) = 433 MeV, being the average M(L = 1) - M(L = 0) » 500 MeV. 
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D sJ (2860), or in the charmonium spectra: the X(3872), the X(3940), the y(3940), the Z(3940), the 
F(4260), and the Z(4430) among others [6]. It seems nowadays unavoidable to resort to higher order 
Fock space components to tame the bewildering landscape arising with these new findings. Four-quark 
components, either pure or mixed with qq states, constitute a natural explanation for the proliferation of 
new meson states [7, 8, 9]. They would also account for the possible existence of exotic mesons as could 
be stable ccnn states, the topic for discussion since the early 1980s [10, 11]. 

All these scenarios suggest the study of qqqq structures and their possible mixing with the qq systems 
to understand the role played by multiquarks in the hadron spectra. The manuscript is organized as 
follows. In Section 2. the variational formalism necessary to evaluate four-quark states is discussed 
in detail with special emphasis on the symmetry properties. In Section 3. the way to exploit discrete 
symmetries to determine the four-quark decay threshold is discussed. In Section 4. the formalism to 
evaluate four-quark state probabilities is sketched. In Section 5. we discuss some examples of four-quark 
states calculated using this formalism. Finally, we summarize in Section 6. our conclusions. 

2. Four-quark spectra 

2.1. Solving the four-body system 

The four-quark {qqqq) problem will be addressed by means of the variational method, specially suited 
for studying low-lying states. The nonrelativistic Hamiltonian will be given by 

H = E + &) + E v ft) , (i) 

i=l ^ i<j=l 

where the potential V(fij) corresponds to an arbitrary two-body interaction. The extension of this 
formalism to consider many-body interactions is discussed in [12, 13]. 

The variational wave function must include all possible flavor- spin-color channels contributing to a 
given configuration. For each channel s, the wave function will be the tensor product of a color (|C Sl )), 
spin (l^)), flavor (\F S3 )), and radial (\R S4 )) component, 

\<t>s) = \C S1 ) <8> \S S2 ) ® |F S3 ) ® \R S4 ) , (2) 

where s = {s±, s 2 , s 3 , s 4 }. The procedure to construct the wave function will be detailed later on. Once 
the spin, color and flavor parts are integrated out the coefficients of the radial wave function are obtained 
by solving the system of linear equations 

E E 0ff K^? I H \ R S) - E I <M = o v j , (3) 

s' s i 

where the eigenvalues are obtained by a minimization procedure. 
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Figure 2. Tetraquark Jacobi coordinates, see Equations (4) for definitions. 




2.2. Four-body wave function 

For the description of the gig 2 g 3 g4 wave function we consider the four-body Jacobi coordinates 
depicted in Figure 2: 



x 



R 



ri - r 2 

r 3 - f 4 
miri + m 2 r 2 
mi + m 2 



(4) 



m 3 r 3 + m4r*4 
m 3 + m 4 



where indices 1 and 2 will stand for quarks and 3 and 4 for antiquarks. Let us now describe each 
component of the variational wave function (2) separately. The total wave function should have 
well-defined permutation properties under the exchange of identical particles: quarks or antiquarks. 
The Pauli principle must be satisfied for each subsystem of identical particles 2 . This imposes restrictions 
on the quantum numbers of the basis states. 

2.3. Color space 

There are three different ways of coupling two quarks and two antiquarks into a colorless state: 



[(31*0 («)] = {|3 12 334),|6 12 6 3 4)} = {|33)J 2 , |66)i 2 } 
[(giftOta)] = {|li 3 l 2 4),|8 13 8 2 4>} = {|ll) c ,|88) c } 
[(gi&XfeftO] = {|li4l 2 3),|8 14 8 23 )} = {|l , l / ) c ,|8 / 8 / ) c } 



(5a) 
(5b) 
(5c) 



2 One should have in mind that if flavor SU (3) symmetry is assumed, u, d, and s quarks are identical particles. 
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being the three of them orthonormal basis. Each coupling scheme allows to define a color basis 
where the four-quark problem can be solved. Only two of these states have well defined permutation 
properties: |33)J 2 , is antisymmetric under the exchange of both quarks and antiquarks, (AA), and 1 66) ,J: 2 
is symmetric, (SS). Therefore, the basis Equation (5a) is the most suitable one to deal with the Pauli 
principle. The other two, Equations (5b) and (5c), are hybrid bases containing singlet- singlet (physical) 
and octet-octet (hidden-color) vectors. The three basis are related through [15, 16]: 

|H)c = y||33)^+y||66)f (6) 

|88) c = -y||33)J 2 + V /||66)f, 

and 

|l'l')c = -\fl\^)l 2 + \fl\M)l 2 (V) 
|8'8') c = yf|33)^ + V /?|66)-. 



Table 1. Color matrix elements. 



6 (Ai-A 2 ) (A 3 -A 4 ) (Xi-A 3 ) (A 2 -A 4 ) (Ai ■ X 4 ) (A 2 • A 3 ) 

J 2 (33|6|33)J 2 -8/3 -8/3 -4/3 -4/3 -4/3 -4/3 

J 2 (66|6|66)J 2 4/3 4/3 -10/3 -10/3 -10/3 -10/3 

i 2 (33|0|66)J 2 -2y/2 -2y/2 2y/2 2y/2 



To evaluate color matrix elements the two-body color operators are introduced in the same manner as 
in angular momentum theory, 

VA, = ^(a 2 .-A 2 -A 2 ), (8) 

where Aj are the SU (3) c Gell-Mann matrices acting on quark i, and A 2 ^ is the Casimir operator. For an 
irreducible representation tp(Xjj?j, the eigenvalue of the Casimir operator is given by: 

A 2 ^(A/i) = ^ (A 2 + /1 2 + X/i + 3A + 3/i) xjj(Xfi) . (9) 
In the color space a quark is described by 3 C = (10) and an antiquark by 3 C = (01), so 

A 2 ^(10) = A 2 [3 c ] = y[3 c ] = ^(10), (10) 

-> 16_ 16 

A^(01) = A 2 [3 c ] = y[3 c ] = y^(01). 

Two quarks in a symmetric state, 6 or 6, have (A/i) = (20) and therefore 

- - 40 

A 2 ^(20) = A 2 [6 C ] = A 2 [6 C ] = T V(20), (11) 
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while two quarks in an antisymmetric state, 3 or 3, have (Xfx) = (01), being the same value as 
Equation (10). Using these expressions, the color matrix elements summarized in Table 1, may be 
easily evaluated. 

2.4. Spin space 

The spin part of the wave function can be written as 

[(siS 2 )5 12 3 S4)s34]s = l 5 12^34)s 2 (12) 

where the spin of the two quarks (antiquarks) is coupled to Si 2 (S 34: ). Two identical spin-1/2 fermions in 
a S = state are antisymmetric (A) under permutations while those coupled to S — 1 are symmetric (S). 
In Table 2 we have included the corresponding vectors for each total spin together with their symmetry 
properties. 



Table 2. Spin basis vectors for all possible total spin states (S). The "Symmetry" column 
stands for the symmetry properties of the pair of quarks and antiquarks. 



s 


Vector 


Symmetry 





\oo)l 2 


AA 




|H>i 2 


SS 






AS 


1 




SA 




in>i 2 


SS 


2 


n)^ 2 


SS 



Using this notation is straightforward to evaluate the four-body spin matrix elements, 



12, 



(SwSul&i ' &j\S[ 2 S' u ) s — 2Sij(Sij + 1) — 3 



Si2,S'Ss 3 4,,S'Js,S' , 



(13) 



for (ij) = (12) or (34) and where c?j is the spin operator acting over particle i. To calculate the other 
spin operators we should reorder the spin wave function [14] 



ySlS2)s 12 {S3S 4 )s 34 



J2(-l) 2Sl2+S2+2s3+S4+l+S V2TTlV2l^^2S 12 + W2S U + 1 



k.i 



S12 s 3 k 
S4 S S34 



52 Si S12 

53 k I 



yS 1 S 3 )lS 2 



s 4 



Now one can calculate the matrix element for the case S\ — s 2 — S3 — S4 — 



12 

s 



' )l 2 



X 



(<S'l2<S'34|o i l ' <?3 I Sl 2 *S34) 

V2S12 + ly/2S' 12 + 1 V2S34 + W 25 ^ + 1 J2( 2k + X )( 2/ + X ) i 2l ( l + 1 

k,l 

1/2 1/2 
1/2 k I 



S 12 1/2 k 
1/2 S S 3A 



S[ 2 1/2 k 
1/2 S S' u 



3] 



1/2 1/2 S[ 2 
1/2 k I 



(14) 



(15) 
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The same can be done for the other spin operators, (B\ ■ c? 4 ), (B 2 ■ 04) and (B 2 ■ B 3 ), using the expressions 
given above. The results are resumed in Table 3. 



Table 3. Spin matrix elements. 



s 




(B\ • B 2 ) 


(B3 ■ Ba) 


(Bi ■ B3) 


(Bo ■ Ba) 


(ai • Ba) 


(Bo ■ as) 




12 <oo 





00)? 


-3 


-3 



















6 




1 


1 


-2 


-2 


-2 


-2 




I 2 (00 





11), 12 








-V3 




v/3 


v/3 




12 <oi 





oi) 12 


-3 


1 
















12 <io 


6 


10)] 2 


1 


-3 














1 


12 <11 


6 


ll) 12 


1 


1 


-1 


-1 


-1 


-1 




12 <oi 


6 


io)^ 2 








1 


1 


-1 


-1 




12 <io 


6 


ll) 12 








V2 


-V2 


-V2 


V2 




12 <oi 


6 


ll), 12 








-V2 


V2 


-V2 


V2 


2 


12 <11 





ll) 12 


1 


1 


1 


1 


1 


1 



Table 4. Pauli-based classification of four-quark states. / indicates that the quark/antiquark 
pair requires the application of the Pauli principle, being the notation (pair of quarks, pair of 
antiquarks). The third and fourth columns contain the recoupling corresponding to bases (5b) 
and (5c). 



(12)(34) 


Pauli 


(13)(24) 


(14)(23) 


(nn)(nn) 


(S,S) 


(nn)(nn) 


(nn)(nn) 


(nn)(nQ) 




(nn)(nQ) 


(nQ)(nn) 


{nn)(Q 1 Q 2 ) 


(/,/ifQi =Q 2 ) 


{nQi)(nQ 2 ) 


{nQ 2 )(nQi) 


{nQi)(nQ 2 ) 


(X,X) 


{nn)(QiQ 2 ) 


(nQ 2 )(Qin) 


(nQ 1 )(Q 2 Q 3 ) 


(X, / if Q 2 = Q 3 ) 


{nQ 2 ){QiQ 3 ) 


(nQ 3 )(Q 1 Q 2 ) 


(QiQ 2 )(Q 3 Q 4 ) 


(/ifg 1 = g 2 ,/ifg 3 = Q 4 ) 


(QiQ 3 )(q 2 q 4 ) 


(QiQ 4 )(Q 2 Q 3 ) 



2.5. Flavor space 

Before discussing the flavor part of the wave function one must specify the required flavor symmetry, 
SU(2) or 577(3). In the former case, u and d quarks are identical whether in the latter, u, d, and s are 
indistinguishable. In the following, n will stand for light u and d quarks and Q for heavy ones, c or b. s 
quarks will be considered heavy if flavor SU (2) is assumed and light otherwise. 

For the flavor part one finds several different possible four-quark states depending on the number of 
light quarks. They can be classified depending on whether they are made of undistinguishable quarks in 
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one of the pairs (and therefore the Pauli principle must be imposed) or not. In following subsections we 
will discuss the important role played by the Pauli principle in the description of the four-quark states 
properties. This classification is illustrated in Table 4. Symmetry properties of the flavor wave function 
are summarized in Table 5. 

Table 5. Symmetry properties of the flavor wave function under the exchange of quarks 
(the same holds for antiquarks). t If flavor £77(3) is assumed, symmetric and antisymmetric 
flavor wave functions with 1 = 1/2 can be constructed, i.e., (us ± su)/V2). 



Flavor 


Symmetry 


nn I = 


A 


nn 1 = 1 


S 


nnl = 1/2+ 


S/A 


QQI = 


S 



The flavor SU(2) matrix elements can be evaluated by means of the same relations shown in 
Section 2.4.. For those corresponding to flavor £77(3) the procedure will require the explicit construction 
of the flavor wave function by means of the SU(3) isoscalar factors given in [15, 16] 3 . As an example 
we evaluate some of the flavor matrix elements needed for the description of heavy-light tetraquarks. 
They can be obtained using the matrix expression of A a , 





( ° 


1 


0\ 




( '° 


—% 


o ^ 








\ 


A x = 


1 








A 2 = 




i 















-1 











°/ 























J 




fo 





1 \ 









-i \ 






( o 


0^ 




A 4 = 











A 5 = 















A 6 = 





1 











°/ 




V 


i 





J 








1 0) 






( 








\ 




I 


1 

^3 





^ 










A 7 = 




1° 




% 


—i 



A 8 = 

/ 




V 






1 






5 









where, following the same convention, quarks and antiquarks are given by, 

u =u= (1,0,0) (17) 
d = J =(0,1,0) 
s =s = (0,0,1). 

The tetraquark flavor wave function corresponding to two light quarks coupled to total isospin I with 
I z = and two heavy antiquarks can be written as 

|V>) = -=[ud+ {-l) I+1 du][ss\ . (18) 
v2 

3 Note there is no universal agreement in the phase convention regarding the isoscalar factor, so mixing different tables 
from different authors should be done with care. 
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A typical flavor operator is 



(19) 



a=l 



where A" are the SU (3) flavor matrices defined above and n axe the isospin Pauli matrices, both acting 
on quark i. So the same expression obtained for the spin operators holds here: 



1 = _> -3 
7=1 -> 1 



Alternatively one can write the flavor matrix element as 



o=l 



wd + (-l) /+1 dii| ^ 

— I 2^ A i A 2 — 

a=l 



y/2 



dw 



a=l L 



du 



v/2 
du 

ud 



A?A5 ud) + (du A"A£ 
A?A£ iid 

^ <j (u|A» (d|A a |d) + (u|A a |d> 
1 + 2(-l) /+1 = 



du 



0=1 



7 = -3 
7 = 1 1 



Other matrix elements of interest are, 



>|A«A^> 



3 

mm = ~ 



2.(5. Radial space 



(20) 



(21) 



(22) 



The most general radial wave function with orbital angular momentum L = may depend on the 
six scalar quantities that can be constructed with the Jacobi coordinates of the system, they are: x 2 , y 2 , 
z 2 , x ■ y, x ■ z and y ■ z. We define the variational spatial wave function as a linear combination of 
generalized Gaussians, 

n n 

\r S4 ) = =TjSrl (23) 



i=l 



where n is the number of Gaussians we use for each color- spin- flavor component. R l depends on six 
variational parameters, a\, b\, c\, d\, e\, and /], one for each scalar quantity. Therefore, any tetraquark 
will depend on6xnxn s variational parameters (where n s is the number of different channels allowed by 
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the Pauli Principle). Equation (23) should have well defined permutation symmetry under the exchange 
of both quarks and antiquarks, 

P 12 (x^-x)< = P X R\ A (24) 

PM(y^-y)K = P y ir S4 , 

where P x and P y are —1 for antisymmetric states, (A), and +1 for symmetric ones, (S). One can build 
the following radial combinations, (P x P y ) = (SS), (SA), (AS) and (AA): 

(25) 



(SS) = 


> R\ 




Exp 


'-a\x 2 - 


dy 2 


-dz 2 


— d l xv — d xz — f l vz) 






+ 


Exp 


\-a\x 2 - 


b l y 2 


-dz 2 


+ d 1 xv — d xz + f l vz) 






+ 


Exp 


'-aix 2 - 


Ky 2 




+ 4xy + d s xz - fiyz) 






+ 


Exp 


'-a\x 2 - 


Ky 2 




- dixy + d s xz + fiyz) 


(SA) = 


> Ri 




Exp 


{-a l s x 2 - 


Uy 2 


-dz 2 


— d 1 xv — d xz — f l vz) 








Exp 


'-a\x 2 - 


Vy 2 


-dz 2 


+ d l xv — d xz + f % vz) 






4- 


1 -/.VI ) 


'-tfx 2 - 


Uv 2 


-dz 2 


4- d l TV 4- d 1 T7 — f % V7\ 








Exp 


{-aix 2 - 


Ky 2 


- clz 2 


- dixy + dgxz + fiyz) 


(AS) = 


>R\ 




Exp 




Ky 2 




- d l s xy- dgxz- fiyz) 






+ 


Exp 


{-aix 2 - 


Ky 2 


~ clz 2 


+ d l s xy-dgxz + flyz) 








Exp 


{-aix 2 - 


Ky 2 


- clz 2 


+ d l s xy + dgxz-flyz) 








Exp 


{-aix 2 - 


Ky 2 


- clz 2 


- dixy + dgxz + f l s y z) 


(AA) = 


>R\ 




Exp 




Ky 2 




- d l s xy- dgxz- f s yz) 








Exp 


{-aix 2 - 


Ky 2 




+ dixy - dgxz + fiyz) 








Exp 


{-aix 2 - 


Ky 2 




+ dixy + dgxz - fiyz) 






+ 


Exp 


{-aix 2 - 


Ky 2 




- dixy + dgxz + f s yz) 



By defining the function 

g(si, -52, s 3 ) = Exp (-aix 2 - bly 2 - d s z 2 - s^xy- s 2 d s xz- s 3 f l s yz) 
we can build the vectors 

/ #(+,+,+) \ 

<?(-+,-) 
<?(--+) 
V <?(+,- -) 



and 



ass 

CiSA 
dAS 
O-AA 



( + ,+,+,+) 
(+,-+,-) 

(+,+,--) 
( + ,-,-,+), 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 
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what allows to write in a compact way Equations (25), (26), (27), and (28), 

(32) 



(SS) = 


> R\ 


= Oiss ■ 


G i 


(SA) = 


> r\ 


= as a 


■ G i 


(AS) = 




= oils 


• G i 

s 


(AA) = 




= ®2a 


■ G i 

s 



Such a radial wave function includes all possible relative orbital angular momenta coupled to L = 0. 
This can be seen through the relation: 



oo oo oo 



Exp {-dixy- e\xz- f s yz) = ^ E E E [[ Y ^)Ye y (y)k Y eM] (33) 



(24 + i)(24 + 1)(24 + l) (4040|4) (4040I4) (40£ 3 o|4) 



£1/2/3 



'i <-z 

h 4 4 



(v^W(^)) (^^WOM 

where 4> 4 an d 4 are the orbital angular momenta associated to coordinates x, y and z, and I a (x) are 
the modified Bessel functions. 

The radial wave functions defined above have also well-defined symmetry properties on the z 
coordinate. Being P(i2)(34)(5* — > — z)-Rg 4 = ^-^s 4 one obtains, 

(34) 



P{12){U)R\ 


= +R\ 


P{12){U)Rt. 


= -R\ 


-P(12)(34)-R 3 


= -R\ 


-P(12) (34)^4 


= +R\ 



To evaluate radial matrix elements we will use the notation introduced in Equation (32): 

(Bt r \f(x,y,z)\R? p )= [ (a s ^Gi)f(x,y,z)(a S 0-G j sl )dV = a s ,-F ij -a 8p , (35) 
Jv 

where 7 and (3 stand for the symmetry of the radial wave function and i™ is a matrix whose element 
(a, b) is defined through, 

Fl = [ {G%{Gi,) b f{x,y,z)dV , (36) 
Jv 

being (G l s ) a the component a of the vector G\. ^From Equation (29) one obtains, 

g(si, s 2 , s 3 )g(s' 1 , s' 2 , 4) = Exp (-a i:j x 2 - b^y 2 - CijZ 2 - s^xy - e i:j xz - f^yz) , (37) 

where we have shortened the previous notation according to and dij = (s\di + 

s[dj). Therefore, all four-body radial matrix elements will contain integrals of the form 

I = Exp (-dijX 2 - b i:j y 2 - Cijz 2 - Sijxy - e^xz - fijyz) f(x, y, z)dxdydz, (38) 
Jv 
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where the functions f(x, y, z) are the potentials. Being all of them radial functions (not depending on 
angular variables) one can solve the previous integral by noting: 



3 ^3 
2 //". \ 2 



jExp[- ^2 A ij x i x j ]f(\^a k x k \)dx 1 ...dx n = y^rjj A7r \^~J F (^ijJ), ( 39 ) 



where 



— - = a- A' 1 - a (40) 

= I e-^/(«)« a d« 

cfetA > 
1 ^ 

One can extract some useful relations for the radial matrix elements using simple symmetry properties. 
Let us rewrite Equation (35) 



Rp xP Pz f(x, y, z)R p , p , p ,dxdydz. 




x J y J z 



If f(x, y, x) depends only in one coordinate, for example x, the integrals over the other coordinates will 
be zero if one of them has different symmetry properties, P y ^ P' or P z ^ P' z in our example. Therefore 

(4|/(:r)|4) oc 5^ (42) 

<i^|/(y)l4) oc 5^ 

<ify/(z)|i#) oc 5 7/3 

(i? 7 | Constant | RV} oc £ 7j g . 

The radial wave function described in this section is adequate to describe not only bound states, but also 
it is flexible enough to describe states of the two-meson continuum within a reasonable accuracy. We 
will came back to this point in Sect. 5. 

2. 7. Parity and C— parity 

The parity of a tetraquark can be calculated as 

—x ' 

(-1)^+4^^ (43) 
or using Equations (24) and (34), 

P [RJJX, y, z)] = P 12 P 3 4^12)(34)<(£, V, Z) = P.PyPz f, i) , (44) 





/ 


x - 


* — X 


\ 


P[<(f,fl5)] = < 




y- 


* -y 






V 


z — 


* — z 


/ 
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what in our case implies 



K x,y,z 



(+)(+)(+)R\ 

R (-)(-) '4 





' +R\ ' 






> = < 


> 






> 





+RUx,y,z) 



(45) 



Hence, this formalism describes positive parity states, being thus adequate to study tetraquark ground 
states. 



Table 6. The action of C over the spin part or the wave function. 



s 


= 


c\oo)l 2 


= +mi 2 






c\n)l 2 


= m)i 2 






c\oi) 12 


= -mi 2 


s 


= 1 


C\10)f 


= -\oi)l 2 








= +lll) s 12 


s 


= 2 


c\n)l 2 


= +|iiL 12 



From Equation (45) one can see that, not only the total wave function will be an eigenstate of the 
parity operator, but also each component will be. This is not the case for C-parity, where only the 
total wave function will be an eigenstate and therefore it must be obtained numerically for each state. 
The tetraquark C— parity will depend on the variational parameters and on the /3s coefficients. This 
dependence is contained in the action of the C— parity operator over the different parts of the wave 
function which will give us the following relations: 

C\Rl 4 (x,y,z))=Rl 4 (y,x,-z), (46) 

and if a\ = b\ and e\ = /j (what is a very common result), 



C\Rl 4 (x,y,z)) 



S4 


= 1 - 


+R\ 


Si 


= 2 - 


--RI 


Si 


= 3 - 


-> —R l 2 


Si 


= 4 - 


+R\ 



(47) 



The action of C over the spin part or the wave function is summarized in Table 6. The action of C over 
the flavor part of the wave function has to be evaluated individually once the wave functions have been 
constructed. 

2.8. qq <-> qqqq mixing 

Many of the possible four-quark systems may present J PC quantum numbers that can be reached not 
only by means of qqqq configurations but also by qq ones and with similar energies. In these cases the 
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possibility of a mixing between them cannot be discarded a priori and therefore, their most general wave 
function will read 

|B = 0) = J^^n \{qq) n ) = fti \qq) + ^2 \qqqq) + •••• (48) 



These particular systems may be described by a hamiltonian 

H = H + H 1 being H = | 



qqqq 



V, 



qq^qqqq 



qq^qqqq 





(49) 



where the nondiagonal terms can be treated perturbatively, therefore allowing to solve the two- and 
four-body sectors separately. 

Figure 3. Coupling between qqqq and qq configurations. 







a 



The Hamiltonian Hi describes the mixing between two- and four-body configurations. Its explicit 
expression would require the knowledge of the operator annihilating a quark-antiquark pair into the 
vacuum. This could be done, for example, using a 3 Po model, but the result will introduce an additional 
degree of uncertainty on the parametrization used to describe the vertex. Such a parametrization is 
determined by the energy scale at which the transition qqqq <-> qq takes place. For the sake of simplicity 
this can be parametrized by looking to the quark pair that it is annihilated, and not to the spectator quarks 
that will form the final qq state: 



V, 



qq^qqqq 



(nnnn\V\nn) = (nsns\V\ss) = (nnns\V\ns) = C n 
(ssss\V\ss) = (nsns\V\nn) = (nsss\V\ns) = C s 



(50) 



A sketch of these mixing interactions is drawn in Figure 3. Such approach has been used in a 
series of papers to describe the light-scalar mesons and the open-charm and open-bottom meson 
sectors [17, 18, 19, 20]. 
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Table 7. Lowest two-meson thresholds in the uncoupled (UN) and coupled (CO) schemes 
for two particular cncn (upper) and ccnn (lower) states, see text for details. They have been 
calculated using the CQC model, see Section 5. for details. M\ M 2 \l indicates the lowest 
threshold and L its relative orbital angular momentum. Energies are in MeV. 





UN 


CO 


s 





1 


2 





1 


2 




I 





1 





1 


J PC 


= 1++ 




J/ipuj\ s 






J/ipp\s 




J/lpu\s,D 


J/lpp\s,D 


(L 


= 0) 




3745 






3838 




3745 


3838 


jPC 


= I" 


DD\ P 


n c u\ P 


D*D*\ P 


DD\ P 


J/ipn\p 


D*D*\ P 


n c u)\ P 


J/tp7l\p 


(L 


= 1) 


3872 


3683 


4002 


3872 


3590 


4002 


3683 


3590 


J p 


= 1+ 




DD*\ S 






DD*\ S 




DD*\ S>D 


DD*\ S>D 


(L 


= 0) 




3937 






3937 




3937 


3937 


J p 


= 1- 


DD\ P 


DD*\ P 


D*D*\ P 


DDxls 


DD*\ P 


D* D*j\ s ,D 


DD\ P 


DD*\ P 


(L 


= 1) 


3872 


3937 


4002 


4426 


3937 


4499 


3872 


3937 



3. Four-quark stability and threshold determination 



The color degree of freedom makes an important difference between four-quark systems and ordinary 
baryons or mesons. For baryons and mesons it is not possible to construct a color singlet using a subset 
of the constituents, thus only qq or qqq states are proper solutions of the two- or three-quark interacting 
hamiltonian and therefore, all solutions correspond to bound states. However, this is not the case for 
four-quark systems. The color rearrangement of Equations (6) and (7) makes that two isolated mesons 
are also a solution of the four-quark hamiltonian. In order to distinguish between four-quark bound 
states and simple pieces of the meson-meson continuum, one has to analyze the two-meson states that 
constitute the threshold for each set of quantum numbers. 

These thresholds must be determined assuming quantum number conservation within exactly the same 
model scheme (same parameters and interactions) used in the four-quark calculation. When dealing with 
strongly interacting particles, the two-meson states should have well defined total angular momentum (J) 
and parity (P), the coupled scheme. If two identical mesons are considered, the spin-statistics theorem 
imposes a properly symmetrized wave function. Moreover, C— parity should be conserved in the final 
two-meson state for those four-quark states with well-defined C— parity. If noncentral forces are not 
considered, orbital angular momentum (L) and total spin (S) are also good quantum numbers, being this 
the uncoupled scheme. 

An important property of four-quark states containing identical quarks, like for instance the QQnn 
system, that is crucial for the possible existence of bound states, is that only one physical threshold 
is allowed, (Qn)(Qn) for the case of heavy-light tetraquarks. Consequently, particular modifications 
of the four-quark interaction, for instance a strong color-dependent attraction in the QQ pair, would 
not be translated into the asymptotically free two-meson state. As discussed in [21], this is not a 
general property of four-quark spectroscopy, since the QQnn four-quark state has two allowed physical 
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thresholds: (QQ)(nn) and (Qn)(nQ). The lowest thresholds for nnQQ states are given in [21], for 
nQnQ states in [22], and those for QQQQ in [23]. We give in Table 7 the lowest threshold for same 
particular cases to illustrate their differences. We show both the coupled (CO) and the uncoupled (UN) 
schemes together with the final state relative orbital angular momentum of the decay products. We would 
like to emphasize that even when only central forces are considered the coupled scheme is the relevant 
one for experimental observations. 

The relevant quantity for analyzing the stability of any four-quark state is A E , the energy difference 
between the mass of the four-quark system and that of the lowest two-meson threshold, 



where E± q stands for the four-quark energy and E(M±, M 2 ) for the energy of the two-meson threshold. 
Thus, < indicates that all fall-apart decays are forbidden, and therefore one has a proper bound 
state. Ae > will indicate that the four-quark solution corresponds to an unbound threshold (two free 
mesons). 

4. Probabilities in four- quark systems 

As discussed in the previous sections four-quark systems present a richer color structure than ordinary 
baryons or mesons. While the color wave function for standard mesons and baryons leads to a single 
vector, working with four-quark states there are different vectors driving to a singlet color state out of 
colorless or colored quark-antiquark two-body components. Thus, dealing with four-quark states an 
important question is whether we are in front of a colorless meson-meson molecule or a compact state, 
i.e., a system with two-body colored components. While the first structure would be natural in the naive 
quark model, the second one would open a new area on the hadron spectroscopy. 

To evaluate the probability of physical channels (singlet-singlet color states) one needs to expand any 
hidden-color vector of the four-quark state color basis in terms of singlet-singlet color vectors. Given a 
general four-quark state this requires to mix terms from two different couplings, Equations (5b) and (5c). 
If (<?i)<?2) or (q 3 , g 4 ) are identical quarks/antiquarks then, a general four-quark wave function can be 
expanded in terms of color singlet- singlet nonorthogonal vectors and therefore the determination of the 
probability of physical channels becomes cumbersome. 

In [24] the two Hermitian operators that are well-defined projectors on the two physical singlet-singlet 
color states were derived, 



A E = E 4q -E(M 1 ,M 2 ), 



(51) 




(52) 



where P, Q, P, and Q are the projectors over the basis vectors (5b) and (5c), 



Q 



P 



ll) c c(ll 

88) c c (88 



(53) 
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and 

P = |l'l') cC (l'l'| 

Q = | 8'8') c c (8'8' | . (54) 

Using them and the formalism of [24], the four-quark nature (unbound, molecular or compact) can be 
explored. Such a formalism can be applied to any four-quark state, however, it becomes much simpler 
when distinguishable quarks are present. This would be, for example, the case of the nQnQ system, 
where the Pauli principle does not apply. In this system the bases (5b) and (5c) are distinguishable due 
to the flavor part, they correspond to [(nc) (en)] and [(nn) (cc)] as indicated in Table 4, and therefore they 
are orthogonal. This makes that the probability of a physical channel can be evaluated in the usual way 
for orthogonal basis [19]. The non-orthogonal bases formalism is required for those cases where the 
Pauli Principle applies either for the quarks or the antiquarks pairs, see Table 4. Relevant expressions 
can be found in [24]. 



Table 8. Mass, in MeV, and flavor dominant component of the light scalar-isoscalar mesons. 



State 


PDG 


Mass 


Flavor 


/o(600) 


400-1200 


568 


(nn 1P ) 


/ (980) 


980±10 


999 


(nnnn) 


/o(1200- 1600) 


1400±200 


1299 


(ss 1P ) 


/o(1370) 


1200-1500 


1406 


(nn 2P ) 


/o(1500) 


1507±5 


1611 


(nsns) 


/o(1710) 


1714±5 


1704 


(glueball) 


/o(1790) 


motf 


1782 


(nn 3P ) 


/o(2020) 


1992±16 


1902 


(ssss) 


/o(2100) 


2103±17 


1946 


(ss 2 p) 


/o(2200) 


2197±17 


2224 


(ss 3 p) 



5. Some selected results 

To illustrate the formalism we have introduced, we discuss some illustrative results. We make use 
of a standard quark potential model, the constituent quark cluster (CQC) model. It was proposed in the 
early 90's in an attempt to obtain a simultaneous description of the nucleon-nucleon interaction and the 
baryon spectra [25]. Later on it was generalized to all flavor sectors giving a reasonable description of 
the meson [26] and baryon spectra [27, 28, 29]. Explicit expressions of the interacting potentials and a 
detailed discussion of the model can be found in [26]. 

The performance of the numerical procedure we have presented described can be checked by 
comparing with other methods in the literature to understand its capability and advantages. Ref. [21] 
makes use of a hyperspherical harmonic (HH) expansion to study heavy-light tetraquarks, obtaining a 
mass of 3860.7 MeV (K max = 24) for the (L, S, I) = (0,1, 0) cenn state using the CQC model. The 
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variational formalism described here gives a value of 3861.4 MeV (with 6 Gaussians), in very good 
agreement. Concerning the unbound states, belonging to the two-meson continuum, the variational 
is able to describe reasonably their energies and root mean square radii. For the unbound (L, S, I) = 
(0, 0, 1) ccnn state the variational method gives a value of A E = +5 MeV to be compared with the value 
obtained with the HH formalism (K = 28), = +33. This is due to the flexibility of the expansion in 
terms of generalized Gaussians and its ability to mimic the oscillatory behavior of the continuum wave 
functions, something that is more difficult using an expansion in terms of Laguerre functions [21]. 

Figure 4. Regge trajectories for the scalar-isoscalar mesons. The squares represent the 
results of Table 8. The lower solid line corresponds to nn systems and the upper line to 
ss systems. The dashed lines correspond to the mass of those states with a large non— qq 
component. 
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Let us now discussed some particular examples where four-quark structures could be present. First 
of all we center our attention on the light scalar-isoscalar mesons. In [18] scalar mesons below 2 
GeV were studied in terms of the mixing of a chiral nonet of tetraquarks with conventional qq states 
using the scheme described in Section 2.8.. We show in Table 8 results for the energies and dominant 
flavor component of the scalar-isoscalar mesons when considering also the mixing with a scalar glueball 
based on intuition from lattice QCD [30, 31, 32, 33]. The results show a nice correspondence between 
theoretical predictions and experiment. This assignment suggests that there are four isoscalar mesons that 
are not dominantly qq states, they are the / (980) (dominantly a nnnn state), the / (1500) (dominantly 
a nsns state), the / (1710) (dominantly a glueball) and the / (2020) (dominantly a ssss state). This 
is clearly seen in Figure 4 where we have constructed the two Regge trajectories associated to the 
isoscalar mesons. As it is observed the masses of the / (600), / (1200 — 1600), / (1370), / (1790), 
/o(2100), / (2200) fit nicely in one of the two Regge trajectories, while those corresponding to the 
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/o(980), / (1500), / (1710), / (2020) do not fit for any integer value. The exception would be the 
/o(2020) that it is the orthogonal state to the / (2100) having almost 50% of four-quark component. 
The glueball component is shared between the three neighboring states: 20 % for the /o(1370), 2 % 
for the / (1500) and 76 % for the / (1710). These results assigning the larger glueball component to 
the /o(1710) are on the line with Refs. [31, 32] and differ from those of Refs. [34, 35, 36] concluding 
that the /o(1710) is dominantly ss and Ref. [37] supporting a low-lying glueball camouflaged within the 
/ (600) peak. 

Table 9. Probabilities (P), in %, of the wave function components and masses (QM), in 
MeV, of the open-charm and open-bottom mesons with 1 = (left) and 1 = 1/2 (right) once 
the mixing between qq and qqqq configurations is considered. Experimental data (Exp.) are 
taken from Ref. [6] . 



/ = 


I = 1/2 


J p = 0+ 


J p = 1+ 


J p = 0+ 


QM 


2339 


2847 


QM 


2421 


2555 


QM 


2241 


2713 


Exp. 


2317.8±0.6 




Exp. 


2459.6±0.6 


2535.4 ±0.6 


Exp. 


2352±50 




P(cnsn) 


28 


55 


P(cnsn) 


25 


~ 1 


Picnnn) 


46 


49 


P(CSl3p) 


71 


25 


P(cSiip) 


74 


~ 1 


P(cn 1P ) 


53 


46 


P(CS 2 3 P ) 


~ 1 


20 


P(CSl3p) 


~ 1 


98 


P(cn 2P ) 


~ 1 


5 


QM 


5679 


6174 


QM 


5713 


5857 


QM 


5615 


6086 


P(bnsn) 


0.30 


0.51 


P(bnsn) 


0.24 


~ 0.01 


P(bnnn) 


0.48 


0.46 


P(bs l3p ) 


0.69 


0.26 


P(6sii P ) 


0.74 


~ 0.01 


P(bn 1P ) 


0.51 


0.47 


P(fe 2 3p) 


~ 0.01 


0.23 


P(6s l3p ) 


~ 0.01 


0.99 


P(bn 2P ) 


~ 0.01 


0.07 



Another interesting scenario where four-quark states may help in the understanding of the 
experimental data is the open-charm meson sector [17, 19, 20]. The positive parity open-charm mesons 
present unexpected properties quite different from those predicted by quark potential models if a pure 
cq configuration is considered. We include in Table 9 some results considering the mixing between 
cq configurations and four-quark states. Let us first analyze the nonstrange sector. The 3 Po cn pair 
and the cnnn have a mass of 2465 MeV and 2505 MeV, respectively. Once the mixing is considered 
one obtains a state at 2241 MeV with 46% of four-quark component and 53% of cn pair. The lowest 
state, representing the L>q(2308), is above the isospin preserving threshold Dn, being broad as observed 
experimentally. The mixed configuration compares much better with the experimental data than the pure 
cn state. The orthogonal state appears higher in energy, at 2713 MeV, with and important four-quark 
component. 

Concerning the strange sector, the D*j(2317) and the _D s j(2460) are dominantly cs J = + and 
J = 1 + states, respectively, with almost 30% of four-quark component. Without being dominant, it is 
fundamental to shift the mass of the unmixed states to the experimental values below the DK and D*K 
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thresholds. Being both states below their isospin-preserving two-meson threshold, the only allowed 
strong decays to D*n would violate isospin and are expected to have small widths. This width has been 
estimated assuming either a qq structure [38, 39], a four-quark state [40] or vector meson dominance [41] 
obtaining in all cases a width of the order of 10 keV. The second isoscalar J p = 1 + state, with an energy 
of 2555 MeV and 98% of cs component, corresponds to the Z) s i(2536). Regarding the £>*j(2317), it has 
been argued that a possible DK molecule would be preferred with respect to an 7 = cnsn tetraquark, 
what would anticipate an I = 1 cnsn partner nearby in mass [42]. The present results support the 
last argument, namely, the vicinity of the isoscalar and isovector tetraquarks. However, the coupling 
between the four-quark state and the cs system, only allowed for the I = four-quark states due to 
isospin conservation, opens the possibility of a mixed nature for the _D*j (23 17), the remaining 1 = 1 
pure tetraquark partner appearing much higher in energy. The I = 1 J = + and J = 1 + four-quark 
states appear above 2700 MeV and cannot be shifted to lower energies. 

Table 10. Heavy-light four-quark state properties for selected quantum numbers. All states 
have positive parity and total orbital angular momentum L = 0. Energies are given in 
MeV. The notation MiM 2 \i stands for mesons Mi and M 2 with a relative orbital angular 
momentum £. P[|33)* 2 (|66);!; 2 )] stands for the probability of the 33(66) components given 
in Equation (5a) and P[| ll) c ( | 88) c )] for the 11(88) components given in Equation (5b). 
Pmm, Pmm*, and Pm*m* have been calculated following the formalism of [24], and they 
represent the probability of finding two-pseudoscalar (Pmm), a pseudoscalar and a vector 
(Pmm*) or two vector (Pm*m*) mesons 



(S,I) 


(0,1) 


(1,1) 


(1,0) 


(1,0) 


(0,0) 


Flavor 


ccnn 


ccnn 


ccnn 


bbnn 


bbnn 


Energy 


3877 


3952 


3861 


10395 


10948 


Threshold 


DD |s 


DD* \ s 


DD* |s 


BB* \ s 


B\B \ P 


A £ 


+5 


+15 


-76 


-111 


-153 


P[|33)J 2 ] 


0.333 


0.333 


0.881 


0.974 


0.981 




0.667 


0.667 


0.119 


0.026 


0.019 




0.556 


0.556 


0.374 


0.342 


0.340 


p[\w) c ] 


0.444 


0.444 


0.626 


0.658 


0.660 


Pmm 


1.000 








0.254 


Pmm* 




1.000 


0.505 


0.531 




Pm*m* 


0.000 


0.000 


0.495 


0.469 


0.746 



We finally tackled an interesting problem in tetraquark spectroscopy, the molecular or compact nature 
of four-quark bound states. This problem requires the determination of probabilities in non-orthogonal 
bases mathematically addressed in [24]. We show in Table 10 some examples of results obtained for 
heavy-light tetraquarks. One can see how independently of their binding energy, all of them present a 
sizable octet-octet component when the wave function is expressed in the (5b) coupling. Let us first 
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of all concentrate on the two unbound states, A E > 0, one with S = and one with S — 1, given in 
Table 10. The octet-octet component of basis (5b) can be expanded in terms of the vectors of basis (5c) as 
explained in the previous section. Then, the probabilities are concentrated into a single physical channel, 
MM or MM* [MM stands for two identical pseudoscalar D (B) mesons and MM* for a pseudoscalar 
D (B) meson together with its corresponding vector excitation, D* (£>*)]. In other words, the octet-octet 
component of the basis (5b) or (5c) is a consequence of having identical quarks and antiquarks. Thus, 
four-quark unbound states are represented by two isolated mesons. This conclusion is strengthened when 
studying the root mean square radii, leading to a picture where the two quarks and the two antiquarks 
are far away, (x 2 ) 1 ^ 2 ^> 1 fm and (y 2 ) 1 ^ 2 ^> 1 fm, whereas the quark-antiquark pairs are located at a 
typical distance for a meson, (z 2 ) 1 / 2 < 1 fm. Let us now turn to the bound states shown in Table 10, 
Ae < 0, one in the charm sector and two in the bottom one. In contrast to the results obtained for 
unbound states, when the octet-octet component of basis (5b) is expanded in terms of the vectors of 
basis (5c), one obtains a picture where the probabilities in all allowed physical channels are relevant. It 
is clear that the bound state must be generated by an interaction that it is not present in the asymptotic 
channel, sequestering probability from a single singlet-singlet color vector from the interaction between 
color octets. Such systems are clear examples of compact four-quark states, in other words, they cannot 
be expressed in terms of a single physical channel. 



Figure 5. Pmm as a function of A E . 
1 -n 1 




-180 -150 -120 -90 -60 -30 
AE (MeV) 

We have studied the dependence of the probability of a physical channel on the binding energy. For 
this purpose we have considered the simplest system from the numerical point of view, the (S, I) = (0,1) 
ccnn state. Unfortunately, this state is unbound for any reasonable set of parameters. Therefore, we bind 
it by multiplying the interaction between the light quarks by a fudge factor. Such a modification does 
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not affect the two-meson threshold while it decreases the mass of the four-quark state. The results are 
illustrated in Figure 5., showing how in the A E — > limit, the four-quark wave function is almost a 
pure single physical channel. Close to this limit one would find what could be defined as molecular 
states. When the probability concentrates into a single physical channel (PmiM 2 — * 1) me system gets 
larger than two isolated mesons [20]. One can identify the subsystems responsible for increasing the size 
of the four-quark state. Quark-quark ((x 2 ) 1 / 2 ) and antiquark-antiquark ((y 2 ) 1 ^ 2 ) distances grow rapidly 
while the quark-antiquark distance ((z 2 ) 1 / 2 ) remains almost constant. This reinforces our previous result, 
pointing to the appearance of two-meson-like structures whenever the binding energy goes to zero. 

6. Summary 

We have presented a detailed analysis of the symmetry properties of a four-quark wave function 
and its solution by means of a variational approach for simple Hamiltonians. The numerical capability 
of the method has been analyzed. We have also emphasized the relevance of a correct analysis of 
the two-meson thresholds when dealing with the stability of four-quark systems. We have discussed 
the potential importance of four-quark structures in several different systems: the light scalar-isoscalar 
mesons and the open-charm mesons. We have also introduced the necessary ingredients to study the 
nature of four-quark bound states, distinguishing between molecular and compact four-quark states. 

Although the present analysis has been performed by means of a particular quark interacting potential, 
the CQC model, the conclusions derived are independent of the quark-quark interaction used. They 
mainly rely on using the same hamiltonian to describe tensors of different order, two and four-quark 
components in the present case. When dealing with a complete basis, any four-quark deeply bound state 
has to be compact. Only slightly bound systems could be considered as molecular. Unbound states 
correspond to a two-meson system. A similar situation would be found in the two baryon system, the 
deuteron could be considered as a molecular-like state with a small percentage of its wave function 
on the AA channel, whereas the H— dibaryon would be a compact six-quark state. When working 
with central forces, the only way of getting a bound system is to have a strong interaction between the 
constituents that are far apart in the asymptotic limit (quarks or antiquarks in the present case). In this 
case the short-range interaction will capture part of the probability of a two-meson threshold to form a 
bound state. This can be reinterpreted as an infinite sum over physical states. This is why the analysis 
performed here is so important before any conclusion can be made concerning the existence of compact 
four-quark states beyond simple molecular structures. 

If the prescription of using the same hamiltonian to describe all tensors in the Fock space is relaxed, 
new scenarios may appear. Among them, the inclusion of many-body forces is particularly relevant. In 
[12, 13] the stability of QQnn and QQnn systems was analyzed in a simple string model considering 
only a multiquark confining interaction given by the minimum of a flip-flop or a butterfly potential in an 
attempt to discern whether confining interactions not factorizable as two-body potentials would influence 
the stability of four-quark states. The ground state of systems made of two quarks and two antiquarks 
of equal masses was found to be below the dissociation threshold. While for the cryptoexotic QQnn 
the binding decreases when increasing the mass ratio mg/m n , for the flavor exotic QQnn the effect 
of mass symmetry breaking is opposite. Others scenarios may emerge if different many-body forces, 
like many-body color interactions [43, 44] or 't Hooft instanton-based three-body interactions [45], are 
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considered. 
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